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Abstract 

In the problem of estimating the lower and upper tail copula we propose two bootstrap procedures 
for approximating the distribution of the corresponding empirical tail copulas. The first method uses a 
multiplier bootstrap of the empirical tail copula process and requires estimation of the partial derivatives 
of the tail copula. The second method avoids this estimation problem and uses multipliers in the two- 
dimensional empirical distribution function and in the estimates of the marginal distributions. For 
both multiplier bootstrap procedures we prove consistency. 

For these investigations we demonstrate that the common assumption of the existence of continuous 
partial derivatives in the the literature on tail copula estimation is so restrictive, such that the tail copula 
corresponding to tail independence is the only tail copula with this property. Moreover, we are able to 
solve this problem and prove weak convergence of the empirical tail copula process under nonrestrictive 
smoothness assumptions which are satisfied for many commonly used models. These results are applied 
in several statistical problems including minimum distance estimation and goodness-of-fit testing. 

Keywords and Phrases: tail copula, stable tail dependence function, multipher bootstrap, minimum 
distance estimation, comparison of tail copulas, goodness-of-fit 
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1 Introduction 



The stable tail dependence function appears naturally in multivariate extreme value theory as a function 
that characterizes extremal dependence: if a bivariate distribution function F lies in the max-domain of 
attraction of an extreme- value distribution G, then the copula of G is completely determined by the stable 
tail dependence function [see e.g. Einmahl et al. (2008)]. The function is closely related to tail copulas 
[see e.g. Schmidt and Stadtmiiller (2006)] and represents the current standard to describe extremal 
dependence [see Embrechts et al. (2003) and Malevergne and Sornette (2004)]. The lower and the upper 
tail copulas are defined by 



Al(x) = lim tG{xi/t,X2/t), A(7(x) = lim t C{xi/t, X2/t), 



(1.1) 



provided that the hmits exist. Here x = {xi,X2) € := [0, oo]^ \ {(00,00)}, C denotes the (unique) 
copula of the two-dimensional continuous distribution function F, which relates F and its marginals -Fi, -F2 
by 



F(x) = C(Fi(xi),F2(x2)) 



:i.2) 



see 



Sklar|([l959j)], and C(u) = ui + ^2 - 1 + C(l - ui, 1 - U2) = P(Fi(Xi) > l-ui,F2{X2) > 1 - U2) 
denotes the survival copula of X = (Xi,X2) ~ F. The stable tail dependence function / and the upper 
tail copula Au are associated by the relationship 

/(x) = + X2 - A{/(x) VxeM+. 

Since its introduction various parametric and nonparametric estimates of the tail copulas and of the stable 
tail dependence function have been proposed in the literature. Several authors assume that the dependence 



function belongs to some parametric family. Coles and Tawn (1994), Tiago de Oliveira (1980) or Einmahl 



et al. (1993) imposed restrictions on the marginal distributions to estimate multivariate extreme value 



distributions. Nonparametric estimates of the stable tail dependence function have been investigated 



in the pioneering thesis of Huang (1992) and by Qi (1997) and Drees and Huang (1998). Schmidt and 



Stadtmiiller (2006) proposed analogous estimates as in Huang (1992) for tail copulas [except for rounding 
deviations due to the fact that F„(F~(x)), with the generalized inverse function F~ , is not exactly equal 
to x] and gave a new proof for the asymptotic behavior of the estimates. More recent work on inference on 



the stable tail dependence function can be found in Einmahl et al. (2008) and Einmahl et al. (2006), who 



investigated moment estimators of tail dependence and weighted approximations of tail copula processes, 
respectively. 

The present paper has two main purposes. First we clarify some curiosities in the literature on tail 
copula estimation, which stem from the fact that most authors assume the existence of continuous partial 



derivatives of the tail copula 


see e.g. 


Huang 


(1992 


), 


Drees and Huang 


(1998) 


(2006 


), 


Einmahl et al. ( 


2006 


), 


de Haan and Ferreira ( 


2006 


), Peng and Qi 



Schmidt and Stadtmiiller 


(2008) or 


de Haan et al. 



(2008) among others]. However, the (lower or upper) tail copula corresponding to (lower or upper) tail 
independence is the only tail copula with this property, because the partial derivatives of a tail copula 
satisfy 

[limt^oo A(l,t) ifx>0 
\0 ifa; = 0, 

where A denotes either A/, or Ajj, see appendix [B] for details. As a consequence we provide a result 
regarding the weak convergence of the empirical tail copula process (and thus also of the empirical stable 



diA{0,x) 



[1.3) 



tail dependence function) under weak smoothness assumptions (see Theorem 2.2 in the following section). 
The smoothness conditions are nonrestrictive in the sense, that in the case where they are not satisfied, 
the candidate limiting process does not have continuous trajectories. 

Note that similar investigations have recently been carried out by Segers ( |2011 ) in the context of nonpara- 
metric copula estimation. In that paper it is demonstrated that many (even most) of the most popular 
copula models do not have continuous partial derivatives on the whole unit square, which has been the 
usual assumption for the asymptotic behavior of the empirical copula process hitherto. Moreover, it is 
shown how the assumptions can be suitably relaxed such that the asymptotics are not influenced. 
The second objective of the paper at hand is devoted to the approximation of the distribution of es- 
timators for the tail copulas by new bootstrap methods. In contrast to the problem of estimation of 
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the stable dependence function and tail copulas, this problem has found much less attention in the lit- 



erature. Recently, Peng and Qi (2008) considered the tail empirical distribution function and showed 



the consistency of the bootstrap based on resampling (again under the assumption of continuous partial 
derivatives). These results were used to construct confidence bands for the tail dependence function. 
While these authors considered the naive bootstrap, the present paper is devoted to multiplier bootstrap 
procedures for tail copula estimation. On the one hand, our research is motivated by the observation that 



the parametric bootstrap, which is commonly applied in goodness-of-fit testing problems [see de Haan 



et al. (2008)], has very high computational costs, because it heavily relies on random number generation 



and estimation [see also Kojadinovic and Yan (2010) and Kojadinovic et al. (2010) for a more detailed 



discussion of the computational efficiency of the multiplier bootstrap]. On the other hand, it was pointed 



out by Biicher (2011) and Biicher and Dette (2010) in the context of nonparametric copula estimation 



that some multiplier bootstrap procedures lead to more reliable approximations than the bootstrap based 
on resampling. 

In Section [2] we briefly review the nonparametric estimates of the tail copula and discuss their main 
properties. In particular we establish weak convergence of the empirical tail copula process under nonre- 
strictive smoothness assumptions, which are satisfied for many commonly used models. In Section [3] we 
introduce the multiplier bootstrap for the empirical tail copula and prove its consistency. In particular, 
we discuss two ways of approximating the distribution of the empirical tail copula by a multiplier boot- 
strap. Our first method is called partial derivatives multiplier bootstrap and uses the structure of the limit 
distribution of the empirical tail copula process. As a consequence, this approach requires the estimation 
of the partial derivatives of the tail copula. The second method, which avoids this problem, is called direct 
multiplier bootstrap and uses multipliers in the two-dimensional empirical distribution function and in the 
estimates of the marginal distributions. Finally, in Section [4] we discuss several statistical applications of 
the multiplier bootstrap. In particular, we investigate the problem of testing for equality between two tail 
copulas and we discuss the bootstrap approximations in the context of testing parametric assumptions 
for the tail copula. Finally, the proofs and some of the technical details are deferred to an appendix. 



2 Empirical tail copulas 

Let Xi , . . . , X„ denote independent identically distributed random variables with distribution function 
F and denote the empirical distribution functions of F and its marginals Fi and F2 by Fn{x) = 
^"^-^ I{Xi < x}, Fni{xi) = Fn{xi,oo) and F„2(x) = F„(oo,X2), respectively. Analogously, we 
define the joint empirical survival function by -Fn(x) = XlILi ^{-^i > ^} the marginal empirical 
survival functions as Fni = 1 — Fni and Fn2 = 1 — Fn2- Following Schmidt and Stadtmiiller (2006) we 
consider the estimators 



\ f \ n fkxi kx2\ , , n- 



/ kxi kx'2 
\ n 



n 



for the lower and upper tail copula, respectively, where fc — )• 00 such that k 
denotes the empirical copula (resp. empirical survival copula), that is 



(2.11 



o(n), and Cn (resp. C„ 



C„(u) = Fn{F-^{ui),F-^{u2)), CnH = F„ (F^l (m ) , (^^2)) 
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Here, G and G denote the (left-continuous) generalized inverse functions of some real distribution 
function G and its corresponding survival function G = 1 — G defined by 



G-{p) :-- 



(mi{x eR\G{x) >p}, 0<p<l 
\sup{x G M| G(x) = 0}, p = 0, 



G-{p) :-- 



sup{x £ M I G{x) >p}, 0<p<l 
mi{x eR\G{x) = 0}, p = 0. 



It is easy to see that the estimators A^, and Ajj are asymptotically equivalent to the estimates 
1 " 

-Y,I{R{Xii) < kxi,R{X^2) < A;x2} = Al(x) + 0(1/A;), 



1=1 



1 " 

-^I{i?(Xa) >n-kxi,R{Xi2) > n - fcxs} = At/(x) + 0(l/fc) 



(2.2) 



(2.3) 



2 = 1 



where R{Xij) = nFni{Xji) denotes the rank of Xij among Xij, . . . Xnj {j = 1, 2), see also [Huang (1992) 



for an alternative asymptotic equivalent estimator. Therefore we introduce analogs of (2.2) and (2.3 
where the marginals Fi and F2 are assumed to be known, that is 



= IEH^i(^^i) ^ — >^2(^.2) < — } 

k ^-^ n n 



i=l 
n 



k ^-^ n n 

i=l 



(2.4) 
(2.5) 



For the sake of brevity we restrict our investigations to the case of lower tail copulas and we assume 
that this function is non-zero in a single point x G [0, 00)^ [and as a consequence non-zero everywhere on 
[0,00)^, see Theorem 1 in Schmidt and Stadtmiiller (2006)]. 

Let Boo (Ml) denote the space of all functions / : — )■ M, which are locally uniformly bounded on every 
compact subset of Ml , equipped with the metric 



d(/i,/2) =£2-^(11/1 -Mb A 1), 



where the sets Tj are defined recursively by Tsi = T^i-i U [0,i]^,r3 



i-l 



732-2 U ([0,i] X {00}),r3i_2 



U ({00} X [0,i]), To = and where ||/||Ti = sup^gTi 1/(^)1 denotes the sup-norm on Tj. Note that 



J 3(2-1) 

with this metric the set ;Soo(M+) is a complete metric space and that a sequence fn in Bo 



converges 



with respect to d if and only if it converges uniformly on every Tj, see Van der Vaart and Wellner (1996). 

p 

Throughout this paper /°°(r) denotes the set of uniformly bounded functions on a set T, — >■ denotes 
convergence in (outer) probability and ^ denotes weak convergence in the sense of Hoffmann-j0rgensen, 
see e.g. 



Van der Vaart and Wellner (1996). 



Schmidt and Stadtmiiller (2006 ) assumed that the lower tail copula A^ satisfies the second-order condition 

Al(x) -fC(xi/t,X2/t) 



lim 

t—>-co 



A{t) 



5(x) 



(2.6) 



locally uniformly for x = (xi, X2) G M^, where 5 is a non-constant function and the function A : [0, 00) — )• 
[0,00) satisfies iiuit^oo A{t) = 0. Under this and the additional assumptions A^, 7^ 0, VkA{n/k) — )• 0, 
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k = k{n) — )• oo, k = o{n), they showed that the lower tail copula process with known marginals defined 
in (2.4) converges weakly in i3oo(]R+), that is 



Vk (Al(x) 



Al(x) ) (b.^^(x) 



(2.7) 



where G^^ is a centered Gaussian field with covariance structure given by 

EG^Jx)G^Jy) = Al(xi Ayi,X2 Aya). 
For the empirical tail copula Ai(x) they established the weak convergence 



a„(x) = Vk (Al(x) - Al(x) 



(x) 



(2.9) 



in BooO^+), provided that the tail copula has continuous partial derivatives. Here the limiting process 



uj^ has the representation 



X 



(x) - 9iAl(x)G^^(xi,oo) - 92Al(x) 



(00,X2 



(2.10) 



The assumption of continuous partial derivatives is made in the whole literature on estimation of stable 



tail dependence functions and tail copulas. However, as demonstrated in (1.3) there does not exist any 



tail copula Al ^ with continuous partial derivatives at the origin (0,0). With our first result we will 
fill this gap and prove weak convergence of the empirical tail copula process under suitable weakened 



smoothness assumptions. For this purpose we will use a similar approach as in Schmidt and Stadtmiiller 



(2006) since this turns out to be also useful for a proof of consistency of the multiplier bootstrap. First 



we consider the case of known marginals. Due to the second order condition (2.6) the proof of (2.7) can 
be given by showing weak convergence of the centered statistic 

/r , , n 



a. 



„(x) := Vk (^Al(x) - '-^C{xik/n,X2k/n] 



2.11 



Lemma 2.1. If Al ^ and the second order condition (2.6) holds with ^/kA{n/k) — t- 0, where k 
k{n) — )■ cxD and k = o{n), then we have, as n tends to infinity 

n 



a„(x) 



'k (Ai(x) — —C(xik/n,X2k/n) 
k 



in ;Boo(IR+)- Here Gi is a tight centered Gaussian field concentrated on Cp{i 



X 2.12 



with covariance structure 



given in (2.8) and p is a pseudometric on the space defined by 

1/2 



p(x,y) 



E 



(AL(x)-2AL(xAy) + AL(x)) 



1/2 



X = (xi,X2), y = (yi,y2), x a y = (xi A yi,X2 A 2/2) and Cp(M^) C ;Boo(K+) denotes the subset of all 
functions that are uniformly p-continuous on every Ti. 



This assertion is proved in Theorem 4 of Schmidt and Stadtmiiller (2006) by showing convergence of 



the finite dimensional distributions and tightness. For an alternative proof based on Donsker classes see 



Remark A. 2 in the appendix. For a proof of a corresponding result for the empirical tail copula process 



with estimated marginals as defined in (2.9) we will use the functional delta method in (2.7) with some 
suitable functional. 
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Theorem 2.2. Let Al ^ be a lower tail copula whose first order partial derivatives satisfy the condition 

dph-L exists and is continuous on {x G | < Xp < oo} (2-13) 
for p = 1,2. // additionally the assumptions of Lemma\2.1\ are satisfied then we have 



a„(x) = fAL(x) - Al(x) 



(x) 



in ;Boo(K+), where the process 
{xGl2_|xpG{0,oo}}. 



is defined in (2.10) and dpAi, p = 1,2 is defined as on the set 



Theorem 2.2 has been proved by Schmidt and Stadtmiiller (2006) [see Theorem 6, therein] under the 
additional assumption that the tail copula has continuous partial derivatives. As pointed out in the 
previous paragraphs there does not exist any tail copula A^ 7^ with this property. 



3 Multiplier bootstrap approximation 

3.1 Asymptotic theory 

In this section we will construct multiplier bootstrap approximations of the Gaussian limit distributions 
G^^ and Gj^^ specified in (2.7) and (2.9), respectively. To this end let be independent identically 
distributed positive random variables, independent of the Xj, with mean /x in (0,oo) and finite variance 
r^, which additionally satisfy ||'^||2,i := Jq° \/ > ^) dx < 00. We will first deal with the case of known 
marginals and define a multiplier bootstrap analogue of (2.4) by 



7i 

where ^„ = Yli=i denotes the mean of ^1, . . . , We have 

«n (X) = E (I - 1) fnAU.) = ^VkCAi - A,), 

where the function fn,x{Ui) is defined by 



(3.1) 



/n,x(Ui) 



{Uii < kxi/n, Ui2 < kx2/n} , 



(3.2) 



(3.3) 



and 



= (Uii ,Ui2); U^p = Fp{X,p) for p = 1, 2. 
Throughout this paper we use the notation 

Gn^G in D 



(3.4) 



for conditional weak convergence in a metric space (D,(i) in the sense of Kosorok (2008), page 19. To be 
precise, (3.4) holds for some random variables Gn = Gri(Xi, . . . , X„, ^1, . . . ^n); G G O if and only if 



sup |E5/i(G„) -E/i(G)| 4 

heSLi(D) 



(3.5) 
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and 



E^h{GnT - ^nHGn)* ^ for every h G BLi{0), 



(3.6) 



where 



BLi{B) = {f 



< l,|/(/3)-/(7)| <d(/3,7)V7,/3eD} 



denotes the set of all Lipschitz-continuous functions bounded by 1. The subscript ^ in the expectations 
indicates conditional expectation over the weights ^ = (^i, . . . , ^„) given the data and /i(G„)* and h{Gn)* 
denote measurable majorants and minorants with respect to the joint data, including the weights ^. The 



condition (3.5) is motivated by the metrization of weak convergence by the bounded Lipschitz-metric, see 



e.g. Theorem 1.12.4 in Van der Vaart (1998). The following result shows that the process (3.2) provides 



a valid bootstrap approximation of the process defined in (2.11) 



Theorem 3.1. If I^l 7^ and the second order condition (2.6) holds with \/~kA{n/k) 
and k = o(n) we have, as n tends to infinity, 



0,k = k{n) 



00 



k(Ai 



A, 



T 



in the metric space Bo 



Since Theorem 



3.1 



states that we have weak convergence of aj^ to conditional on the data Ui, it pro- 



vides a bootstrap approximation of the empirical tail copula in the case where the marginal distributions 
are known. To be precise, consider S G N independent replications of the random variables ■ ■ ■ ,£,n and 
denote them by .^i^;,, . . . , £,n,b- Compute the statistics a^^, = a^{S,ifi, ■ ■ ■ , 6.n,b) (b = I, ■ ■ ■ , B) and use the 
empirical distribution of a™]^ , • . . , ^ approximation for the limiting distribution of Gj^^ . 
Because in most cases of practical interest there will be no information about the marginals one cannot 



use Theorem 3.1 in many statistical applications. We will now develop two consistent bootstrap approxi- 



mation for the limiting distribution of the process (2.9) which do not require knowledge of the marginals 



Intuitively, it is natural to replace the unknown marginal distributions in (3.1) by their empirical coun- 
terparts, that is 



Ai'-(x) = 
which yields the process 

/3„(x) = '^Vk (Ai'- - Al 



I E hi^^' ^ F-,(kxi/n),X,2 < F-,{kx2/n)} 



(3.7) 



^k 



=1 



in 



1 l{X,i < F-^{kxi/n),Xi2 < F-^{kx2/n)]. 



Unfortunately, this intuitive approach does not yield an approximation for the distribution of the process 



Theorem 3.2. Suppose that the assumptions of Theorem\2.2\ hold. Then we have, as n tends to infinity 



/3n 



Al) ^ G 



in the metric space Boo(!i 
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Although Theorem 3.2 provides a negative result and shows, that the distribution of can not be used 



for approximating the limiting law G^^, it turns out to be essential for our first consistent multiplier 
bootstrap method. To be precise, we note that the distribution of /3„ can be calculated from the data 
without knowing the marginal distributions. As a consequence, we obtain an approximation for the 
unknown distribution of the process Gi .In order to get an approximation of Gi we follow IRemillard 



and Scaillet ( 2009 ) and estimate the derivatives of the tail copula by 



dpAL{x + {h- Xp)ep) ,Xp < h 
,Xp = oo 

where denote the pth unit vector {p = 1,2) and h ~ A:~^/^ tends to with increasing sample size. We 
will show in the Appendix (see the proof of the following Theorem in Appendix A) that these estimates 
are consistent, and consequently we define the process 

<'^™(x) = /3„(x)-a;Ai(x)/3„(:Ei,oo)-5^(x)/3„(oo,X2). (3.8) 

Note that a^"^ only depends on the data and the multipliers ^i, . . . ,^ri- As a consequence, a bootstrap 
sample can easily be generated as described in the previous paragraph and we call this method partial 
derivatives multiplier bootstrap (pdm-hootstrap) in the following discussion. Our next result shows that 
the pdm-bootstrap provides a valid approximation for the distribution of the process G^ . 



Theorem 3.3. Under the assumptions of Theorem 2.2 we have 



in the metric space Soo(K+). 

It turns out that there is an alternative valid multiplier bootstrap procedure in the case of unknown 
marginal distributions, which is attractive because it avoids the problem of estimating the partial deriva- 
tives of the lower tail copula. This method not only introduces multiplier random variables in the two- 
dimensional distribution function but also in the inner estimators of the marginals. To be precise define 

F|(x) = - V ^l{Xii < Xl,Xi2 < X2} 



1 " t 

C«'«(u) = F|(F«-(ni),F4-(n2)). 



and consider the process 



Ai'«(x) := Qx) = 1 < F«-(A:xi/n),X.2 < F«-(fex2/n)} (3.9) 

Throughout this paper we will call this bootstrap method the direct multiplier bootstrap (dm-bootstrap) . 
Theorem 3.4. Under the assumptions of Theorem \2.S\ we have 

atr{^) = -Vk(k^j^{^)-kd^) 



G^^ in Bo 



-))• 



(3.10) 
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3.2 Finite sample results 

In this section we will present a small comparison of the finite sample properties of the two bootstrap 
approximations given in this section. We also study the impact of the choice of the parameter k on the 
properties of the estimates and the bootstrap procedure. For the sake of brevity we only consider data 
generated form the Clayton copula with a coefficient of lower tail dependence = 0.25. The Clayton 
copula is defined by 

C{u;e) = (uf + u^'^ -iy^^\ 9>0, (3.11) 
and is a widely used for modeling of negative tail dependent data. Its lower tail copula is given by 

Ar(x) = 1x7" + x-"' ^'^ 



In Tables [T] and [2] we investigate the accuracy of the bootstrap approximation of the covariances of the 
limiting variable G^^^ . 

We chose three points on the unit circle {e^"^ ^ip = I'Kj'i with I = 1,2,3} and present in the first four 
columns of Table |Tj the true covariances of the limiting process G^^ . The remaining columns show the 
simulated covariances of the process on the basis of 5 • 10^ simulation runs, where the sample size is 
n = 1000 and the parameter k is chosen as 50. This choice is motivated by the left panel of Figure [T] where 
we plot the sum of the squared bias, the variance and the mean squared error of the estimators Ai(e*^'^/'^) 
for AL{e'^P^/^) {£ = 1,2,3). The MSB is minimized for values of k in a neighbourhood of the point 50. 
Note also that the literature provides several data-adaptive proposals for the choice of the parameter k, 



see for example Drees and Kaufmann (1998) oder Gomes and Oliveira (2001) in the univariate context. 



Table [T] now serves as a benchmark for the multiplier bootstrap approximations of the covariances stated 
in Table [2j where we investigate the quality of the approximation by various bootstrap methods. The 
distribution of the multipliers in the dm and pdm bootstrap procedure has been chosen according to 



Biicher and Dette (2010) as = 0) = = 2) = 0.5, such that n = r = 1. For the sake of completeness 



we also investigate the resampling bootstrap considered in Peng and Qi| ( 2008| ) [which is hereafter denoted 



by al^'^]. The estimated covariances given in the first part (lines 3-5) of Table [2] have been calculated by 
1000 simulation runs, where in each run the covariance is estimated on the basis oi B = 500 bootstrap 
replications. The second part (lines 6-8) of Table [2] shows the corresponding mean squared error. 
As one can see all bootstrap procedures yield approximations of comparable quality. Considering only 
the bias in Table [2] the pdm-bootstrap has slight advantages in all cases, while there are basically no 
differences between the dm- and the resampling bootstrap. A comparison of the mean squared error in 
Table [2] shows that the pdm-bootstrap has the best performance on the diagonal. On the other hand, it 
yields a less accurate approximation for the off-diagonal covariances. In this case, the dm-bootstrap yields 
the best results. 

In the right panel of Figure [T] we investigate the sensitivity of the accuracy of the estimators for the 
covariances with respect to the choice of k. For this purpose we calculated the sum of the MSE- values 
given in Table [2] (as well as the variance and squared bias) for various choices of k. As one can see the 
best choices for k lie in an interval of approximate length of 100 around the center k = 200. Compared 
to the "best" value /c = 50 for estimating A^ the optimal values for estimating the covariances of Gj^^ 
are approximately four times larger for both the pdm- and the dm-bootstrap. This increase may be 
explained by the fact that the large bias of A2,(x), A^''(x) and A|^'^(x) for estimating Ai(x) cancels out if 
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True 






TV o TT q TT 

8 ^8 "^8 


8 ^8 "^8 


TT 

8 

2- 
^8 

•^8 


0.0874 0.0754 0.0516 
0.1160 0.0754 
0.0874 


0.0889 0.0737 0.0476 
0.1218 0.0741 
0.0892 



Table 1: Left part: True covariances ofG^^ for the Clayton Copula with Xl = 0.25. Right part: sample 
covariances of the empirical tail copula process an with sample size n = 1000 and parameter A; = 50. 





prfm 
an 


"n 






8 ^8 "^8 


8 ^8 "-'8 


8 ^8 "-'8 


7r 
8 

2- 
^8 

■^8 


0.094 0.072 0.046 
0.130 0.072 
0.094 


0.100 0.071 0.045 
0.136 0.707 
0.099 


0.100 0.070 0.043 
0.136 0.070 
0.099 


7r 
8 

2- 
^8 

■^8 


3.67 4.68 3.65 
8.11 4.87 
3.70 


3.86 3.49 2.72 
8.89 3.25 
3.77 


4.21 3.85 3.21 
8.73 3.64 
3.90 



Table 2: Averaged sample covariances (lines 3-5) and mean squared error x 10^ (lines 6-8) of the 
Bootstrap approximations a?f"^, and ofGx under the conditions of Table 



1 



Clayton (lambda=1/4) 



Clayton (lambda=1/4) 











Bias^2 






- - Var 






MSE 




ft 
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-- Var (dm) 






MSE (dm) 






Bias'^2 (pdm) 






Var (pdm) 






MSE (pdm) 




\\ 




// 
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N -s. . ' ' 

















20 
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Figure 1: Left panel: Averaged MSE, variance and squared bias for the estimation of AL{e^^^^^) (£ = 
l,2,3j by its empirical counterpart AL{e^'^^'^) against the parameter k. Right panel: Averaged MSE, 
variance and squared bias for the bootstrap estimation of the covariances of G i 
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one calculates the difference A^' (x) — Al{'x.) or A^'^(x) — Ai(x). As a result we may choose larger values 
of k which yields a notable decay of the variance. These results indicate that a larger value of k should 
be used in the bootstrap procedure. 

Finally, a comparison of the variance and the bias of the two bootstrap procedures investigated in Figure [T] 
reveals that the pdm bootstrap has a smaller bias than the dm-bootstrap, whereas the variance is slightly 
larger. On the other hand the differences with respect to the MSE are nearly not visible. 



4 Statistical applications 

In this section we investigate several statistical applications of the multiplier bootstrap. In particular we 
discuss the problem of comparing lower tail copulas from different samples, the problem of constructing 
confidence intervals and the problem of testing for a parametric form of the lower tail copula. 



4.1 Testing for equality between two tail copulas 

Let Xi, . . . , X„j and Yi, . . . , Y„2 denote two independent samples of independent identically distributed 
random variables [we will relax the assumption of independence between the samples later on] with 
continuous cumulative distribution function F = C{Fi, F2) and H = D{Hi, H2), respectively. We assume 
that for both distributions the corresponding lower tail copulas, say A^^x and A^^y, exist and do not 
vanish. We are interested in a test for the hypothesis 

no ■■ Al,x = Al,y vs. ni : Al,x / Al,y (4.1) 

Due to the homogeneity of tail copulas we have AL^t'x.) = tAi(x) for all f > 0,x G [0, 00)^, and the 
hypotheses are equivalent to 

^0 : £'(Al,x, AL,y) = vs. Tii : q{Al,x,Al,y) > 0, 
where the distance g is defined by 

q{Al,x,Al,y) ■= {AL,x{cosip, simp) - ALxicosip, simp)) dip (4.2) 

Jo 

^ \aI^{p) - Alyiip))' dip 



and we have used the notation A£j^((^) = AL^x{<^osip,sinp), Af y((/?) = A^^y (cos sin (/?). We propose 
to base the test for the hypothesis (|4.1|) on the distance between the empirical tail copulas and define 



eCAL,x,KY) = / iKxi^) - KYi^)fd^, 



^ kik2 . , kik2 

= 1 — r-r0i^L,x,AL,Y) = , , , 

ki + K2 ki + k2 Jo 

where A£j^((^) = Ai^x(cos((^), sin((/7)), A^ y = A^^y (cos((^), sin((^)) denote the empirical tail copulas 
Aj^^x and A^^y with corresponding parameters ki and k2, satisfying 

kp — 00, kp = o{np) (p = 1, 2) and ki/{ki + k2) — )■ A G (0, 1). 
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We assume that the tail copulas A^^x and A^^y satisfy a second order condition as in (2.6) (with A replaced 
by Ap) and that kp is chosen appropriately, i.e. \/kpAp{kp/np) = o(l). Under the null hypothesis (4.1) of 
equality between the tail copulas we have Sn = Tn with 

rTv/2 

<S^(cos sin ip) d(f, 



where 



Tn 



A;i(Al,x(x) - Al,x(x)) 



k^ 



A:2(AL,y(x) - AL,y(x)). 



ki+k2 ' ' V ^1 + ^2 

Since the two samples X and Y are independent we obtain, independently of the hypotheses, that 



£n \/l - AG: 



£. 



(4.3) 



in the metric space Boo{^\)-, where the stochastically independent two-dimensional centered Gaussian 
fields ^ and ^ are defined in (2.10). This yields by the continuous mapping theorem 

Tn ^ j £ {cos If, sin If) dip =: T 
Jo 

under both the null hypothesis and the alternative. Note that q{^l,x, -^Lx) ~^ q{^l,x, ^l,y), which 
vanishes if and only if the null hypothesis (4.1) is satisfied. Therefore, we can conclude that 



Sn -^-Ho T, Sn —^Hi CO, 



(4.4) 



which shows that a test, which rejects the null hypothesis (4.1) for large values of Tn is consistent. 
In order to determine critical values for the test we approximate the limiting distribution T by the 
multiplier bootstrap proposed in Section [3j For this purpose we exemplarily consider the pdm-bootstrap 
(the extension to the dm-bootstrap is straightforward) using the definition in equation (3.10) and denote 
for any 6 E {1, . . . , B} by ^i^;,, . . . , ^m.fe) Ci,b, • • • i Cn2,b independent and identically distributed non-negative 
random variables with mean /ii (resp. fi2) and variance (resp. r|). We compute for each b and both 
samples the bootstrap statistics as given in (3.8), i.e. 

«X,ni,feW = /3x,ni,6(x) - OiAl^X {^) f^X ,m,bixi, Oo) - ^2 A^^x (x)/3x,ni ,6(oO, X2) , 



a 



pdm 
Y,n2, 



.b(x) = PY,n2,b{^) - 9lAL_y(x)/3y^„2^b(xi,00) - d2ALx{^)PY,n2,b{^^^2), 



where 



^X,ni,fe(x) 

/5y,n2,fe(x) 



fJ-1 1 



E 

i=l 

E 



Ci,b 



1 

T2 \/k2 



i=l 



Ci,b 



ij ^X^l < F„~^i(A:iXi/ni),X,2 < F;;^^{kiX2/ni)}, 
1 ) HYii < H-^^{k2Xi/n2),Y,2 < H-^^{k2X2/n2)], 



and dpA^^x and dpKiy are the corresponding estimates of the partial derivatives {p = 1,2). For all 
X G 1^ and ah 6 G {1, ... , B} define 



^(^'^-'^)(x) :=y^ 



-a 



+ fc2"^ 



pdm 



m / \ 



ki + k2 



-ot 



pdm 
Y,n2,b 



(x), 
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Jo 

By Theorem 3.3 and Theorem 10.8 in Kosorok| (2008), it follows for every b € {1, . . . , B} 



where T"^''-* is an independent copy of T [note that we consider the processes g^'^'^'^^ [n the Banach space 
/°^([0, l]'^)]. From (4.4) we therefore obtain a consistent asymptotic level a test for the null hypothesis 



( |4.1[ ) by rejecting TIq for large values of 5„, that is 

Sn > qf:: (4.5) 

where denotes the (1 — a) quantile of the distribution function Kff'^{s) = Z^^i ^{Tn^'^^''^^ < s}. 
The discussion present so far refers to two independent samples. Nevertheless it is easy to check that 
the methodology of the previous sections also applies if we are faced with paired observations, i.e. Xj is 
not independent of Yj, but ni = n2 = n. In that case we have to set = £,i,b for all i = 1, . . . ,n and 
b = 1, . . . , B. To see this, set Zj = (Xji, Xi2, Yji, Yi2) and denote the (empirical) copula of Zj by (Cn) C. 
Clearly, 

C{ui,U2) = C{ui,U2, 1, 1), D{vi,V2) = C(l, l,Vi,V2), 

C„(ui,U2) = Cn{ui,U2, 1, 1), Dnivi,V2) = C„(l, l,Vi,V2). 

If we set AL,z(x,y) = linii^^oo * C(x/t, y/t), AL,z{x,y) = ^C„(^, ^), we obtain 

Al,x(x) = Al,z(x, oo, oo), Ai^y (y) = Al,z{oo, oo, y), 

Al,x(x) = Al,z(x, oo, oo), AL,y (y) = Al,z(oo, oo, y). 

Under a second order condition on the joint tail copula A^^z the asymptotic properties of the process A^^z 
can be derived along similar lines as given before and the details are omitted for the sake of brevity. As 
the only difference to the preceding discussion note that the occurring limiting fields ^ and G^^^ ^ are 
not independent anymore. Since the asymptotic behavior of the multiplier bootstrap approximations can 
be shown to reflect this dependence we still obtain consistency of the test, the details are again omitted. 
For an investigation of the finite sample property we consider two independent samples of independent 



identically distributed random variables with Clayton copula, see (3.11), with a coefficient of lower tail 
dependence varying in the set {0.25,0.5,0.75}. 



In Tablejsjwe show the simulated rejection probabilities of the pdm- and dm-bootstrap test defined in (4.5 ) 



for various nominal levels on the basis of 1000 simulation runs. The sample size is ni = n2 = n = 1000 
and B = 500 bootstrap replications with ^({0, 2})-multipliers [i.e. P(C = 0) = P(C = 2) = 0.5, such that 
fj, = T = 1] have been used. The parameter k is chosen as either A; = 50 or A; = 200 as suggested by the 
discussion in the preceding paragraph. 

We observe that the nominal level is well approximated by the pdm bootstrap if the coefficient of tail 
dependence is not too large. For a larger coefficient the test tends to be conservative. It is worthwhile 
to mention that the approximation of the nominal level is rather robust with respect to the choice of k. 
A comparison of the performance of the two bootstrap procedures shows that the dm bootstrap test is 
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0.25 
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0.5 
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1 


1 


1 


1 


1 
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0.25 


0.75 


1 
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1 


1 


1 


1 



Table 3: 



Simulated rejection probabilities of the bootstrap tests defined in (4.5) for the hypothesis (4.1) . 



slightly more conservative and this effect is increasing with the coefficient of tail dependence. 
The alternative of different lower tail copulas is detected with reasonable power, where both tests yield 
rather similar results with slight advantages for the pdm-bootstrap. An investigation of the impact of the 
choice of the parameter k under the alternative shows some advantages for k = 200. This may again be 
explained by the fact that bias terms cancel out if one calculates the difference x — 



4.2 Bootstrap approximation of a minimum distance estimate and a computationally 
efficient goodness-of-fit test 

In this section we are interested in estimating the tail copula of X under the additional assumption that 
it is an element of some parametric class, say C = {Al{-;9) \ 9 G ©}. Recently, estimates for parametric 



classes of tail copulas and stable tail dependence functions have been investigated by de Haan et al. (2008) 



and Einmahl et al. (2008) who proposed a censored likelihood and a moment based estimator, respectively. 
In the present section we investigate a further estimate, which is based on the minimum distance method. 
To be precise let Al denote an arbitrary lower tail copula and Al{- -jO) an element in the parametric class 



£ and consider the parameter corresponding to the best approximation by the distance g defined in (4.2) 



0B = T{Al) 



arg min£)(AL, Al( 

See 



(4.6) 



where g is defined in (4.2). We call 6. 



WD 



T(Al) a minimum distance estimator for 9, where A^ is the 
empirical lower tail copula defined in (2.1). Note that 9b is the "true" parameter if the null hypothesis is 
satisfied. 

Throughout this subsection let Xi , . . . , X„ denote independent identically distributed bivariate random 
variables with cumulative distribution function F = C{Fi, F2) and existing lower tail copula A^,. Further- 
more, assume that the standard conditions of minimum distance estimation are satisfied. For a precise 



formulation of these conditions and a proof of the following result we refer to Biicher (2011). 
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Theorem 4.1. If the true tail copula Al satisfies the first order condition (2.13) of Theorem 2.2, then 
the minimum distance estimator is consistent for the parameter Ob corresponding to the best approx- 
imation with respect to the distance g. Moreover, 

G^^ := -eB) = Vkj le,{v) [kii^) - Af(v^)) d^ + op(l) 

^e,{'f)Gl^{ip)d^ 

where Aj^{if) = (cos sine/?), Af = 1^(00899,311199), 70^(92) = ^el^ig^^^' ^ei'f) = ^eA^ (cos 99, sin 93, 61), 
(99) = Gy^^ (cos 99, sin 99) and 



Ai{ip))d^ 



with Aj^{ip;d) = A^, (cos 99, sin 99; 0), The limiting variable Q^^^ is centered normally distributed with vari- 



ance 



a 



[0,7r/2]2 



7es i^hoB {^')r{cos ip, sin 99, cos 99', sin ip') d{ip, ip'), 



where r denotes the covariance functional of the process Gi defined in (2.10) 



In order to make use of the latter result in statistical applications one needs the quantiles of the limiting 
distribution. We propose to use the multiplier bootstrap discussed in the previous section. The following 
theorem shows that the pdm and dm bootstrap yield a valid approximation of the distribution of the 
random variable Q^^^ . 



Theorem 4.2. If the assumptions of the Theorems 3.3 and 3.4 hold and r„ denotes either the process 
^pdm ^y/jg(jj™gy^ ^j^'^ or afK^ (Theorem 3.4) obtained by the pdm- or dm-bootstrap, respectively, then 



r\MD,m 



7^Mz,(^)r^(9^)d9' ^ e^'^, 



where r^{<f) = r„ (cos 99, sin 99), j^md = A-^^^ ^Lioif) (^^^ 

HMn - I 54,,((^)<^4,,(95)^ + a,5|,,,(99)(Af(99;^f«)-Af(99))d9,. 

On the basis of this result it is possible to construct asymptotic confidence regions for the parameter 9 
as well as to test point hypotheses regarding the parameter. In Table |4] we present a small simulation 
study regarding the finite sample coverage probabilities of some confidence intervals for the parameter of 
a Clayton tail copula. This interval is defined as Kli^a = [^^^ — k~^^'^Qi-a/2y^n^^ ~ ^~^^'^Qa/2]^ where 
q/s denotes the estimated /3-quantile of the distribution of ©^'^'^ based on the bootstrap approximation 
provided by Theorem |4.2[ The sample size is n = 1000 and B = 500 bootstrap replications are used for 
the calculation of the quantiles. All coverage probabilities are calculated by 1000 simulation runs. The 
parameter of the Clayton tail copula is chosen such that the tail dependence coefficient varies in the set 
{1/4,2/4,3/4}. 
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k 


50 


200 


50/200 


Al 


90% 95% 


90% 95% 


90% 95% 


0.25 


0.895 0.955 


0.014 0.044 


0.830 0.915 


0.5 


0.893 0.936 


0.779 0.882 


0.888 0.934 


0.75 


0.838 0.887 


0.900 0.949 


0.863 0.894 



Table 4: Simulated coverage probability of the confidence intervals based on the pdm-bootstrap for n = 
1000. In the last two columns the parameter k is chosen as 50 for the estimation of 6, whereas k is chosen 
as 200 for the bootstrap approximation of G s 



In order to investigate the impact of the choice of k we chose three different scenarios: k = 50, k = 200 and 
for the last setting we used two different values of k, namely fc = 50 for the estimator and k = 200 



for the Bootstrap-estimator of the quantiles qf^. This choice is motivated by the findings in Section 3.2 
which indicate that a smaller value of k should be used in the estimator A^,. 

The tables reveal that there is no unique "optimal" choice for k. For = 0.25 the best results are 
obtained for /c = 50 followed by the case of two different values of k [these findings may be compared to 

For k = 200 the large bias of 0^^^ (compare the left-hand side of Figure ll 



3.2 



the results of Section 

entails that the true parameter does not lie in the estimated confidence interval for more than 95% of 
the repetitions. For stronger tail dependence = 0.5,0.75 the choice k = 200 yields better results, with 
almost perfect coverage probabilities for Al = 0.75. It is also of interest to note that the choice /c = 50 
in the estimator A/, and k = 200 in the corresponding bootstrap statistic does not yield an improvement 
with respect to the approximation of the coverage probability compared to the case A; = 50. 

Remark 4.3. As pointed out at the beginning of this section there exist two alternative estimators 



for parametric classes of tail copulas, de Haan et al. (2008) proposed a censored maximum likelihood 



estimator and proved weak convergence to a normal distribution, which involves the partial derivatives 



up to the sixth order of the stable tail dependence function. Einmahl et al. (2008) proposed a method 



of moment type estimator and proved a similar statement as given in Theorem 4.1 for the minimum 
distance estimate. In Table [5] we compare the asymptotic variances of the method of moment and the 
minimum distance estimator for the parameter 6 in the Clayton family chosen such that the coefficient of 
tail dependence A varies in the set {0.1, . . . , 0.9}. The calculated values E\ are defined as 



Asymptotic variance of the minimum distant estimate 
Asymptotic variance of the moment type estimate 



(note that we were not able to obtain the asymptotic variances for the censored maximum likelihood 
estimator, because of the complicated structure of the limiting distribution). The method of moment 



estimator requires the specification of a function g, which was chosen as in [Einmahl et alT (2008) as the 
indicator of the set {x G [0, 1]^ : xi + X2 < 1}. We observe that none of the two estimates is globally 
preferable to the other. For small amounts of tail dependence the minimum distance estimate performs 
slightly better while for increasing tail dependence the moment type estimator is more qualified from an 
asymptotic point of view. 

It is also notable that the dm- and pdm-bootstrap can be used to construct a consistent approximation 



of the asymptotic distribution of the censored likelihood and moment estimator investigated in de Haan 
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A 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


Ex 


0.86 


0.87 


0.89 


0.92 


0.97 


1.03 


1.14 


1.38 


2.11 



Table 5: Relative efficiency of the minimum distant estimate to the moment type estimate. 



et al. 


(2008) and 


Einmahl et al. 


(2008 



distribution of the method of moments and the minimum distance estimator can be represented in the form 
^{Gj;^^, Al, BAl) for some appropriate functional ^ depending on the method of estimation. Here dA^ 
denotes any vector of partial derivatives Al with respect to its coordinates or the parameter. Given the 
functional $ is suitable smooth the bootstrap approximation is then obtained by $(a„, A^,, 9Ai) where 
an is aff"^ of ai^ and OAl is a consistent estimate of OAl- 



In the following we will use the multiplier bootstrap to construct a computationally efficient goodness-of-fit 
test for the hypothesis that the lower tail copula has a specific parametric form, i.e. 



■Ho:ALeC = {ALi-,e)\eGe}, -Hi-.Al^C. 



(4.7) 



This problem has also been discussed in de Haan et al. (2008) and Einmahl et al. (2008) who proposed a 
comparison between a nonparametric and a parametric estimate of the lower tail copula by an L^-distance. 
In both cases the limiting distribution of the corresponding test statistic under the null hypothesis depends 
in a complicated way on the process Gi and the unknown true parameter 9b- While Einmahl et al 



(2008) do not propose any bootstrap approximation, de Haan et al. (2008) proposed to use the parametric 



bootstrap. However, it was pointed out by Kojadinovic and Yan (2010) or Kojadinovic et al. (2010) that 



for copula models, approximations based on multiplier bootstraps are computationally more efficient, 
especially for large sample sizes. We will now illustrate how the multiplier bootstrap can be successfully 



applied in the problem of testing the hypothesis (4.7). 



To be precise, we propose to compare a parametric [using the minimum distance estimate 0^^^] and a 



nonparametric estimate of the tail copula and to reject the null hypothesis (4.7) for large values of the 
statistic 



G0Fn:=kQiAL,ALi-,9^"')) = k I ( A^(^) - A^((^; ^ 



dip, 



where denotes the minimum distan ce estimate. If the standard assumptions of minimum distance es- 

(2011 ) for details] we obtain for the process Hn = Vk (Al - Al{-; 9^^ 



timation are satisfied [see 



Biicher 



under the null hypothesis Hq : Al = Al(-; 9b) 

Hn = Vk (Ai -Al- ^ - 0)) + op(l) 

= Vk (^Al - Al - Se j 7e(^)(Af ((/?) - Af (yp) d^ + 



Under the alternative hypothesis we get an additional summand 
Hn = Vk (Al -Al- <5,(^f ^ -9)- {Al{- 



Al) +op(l). 
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which converges to either plus or minus infinity whenever Al{x,9b) 7^ A/,(x). The continuous mapping 
theorem yields the following result. 



Theorem 4.4. Assume that assumptions of Theorem 4-1 are satisfied. If the null hypothesis is valid then 

2 



dip, 



(4i 



while under the alternative GOFa = J{H^{ip)}^ d(p — )• 00. 

The critical values of the test, which rejects the null hypothesis for large values of GOFn can be calculated 



on the basis of the following theorem. For a proof see Biicher (2011) 



Theorem 4.5. // the assumptions of the Theorems 3.3 and \3.4\ hold and r„ denotes either the process 
^pdm ^ji/jgQ^g^ ^j^^ gj, ^dm ^T^/jg^j-g^ obtained by the pdm- or dm-bootstrap, respectively, then it 



holds independently of the hypotheses that 



Hn ■= r„ - SgMD / IgMoif) r„((/7) dip - G 



MD 



Therefore GOFj;;' = J {Hj;^'^ {p)}'^ dp ^ Z, where Z is defined in d^S]). 



In order to investigate the finite sample properties of a goodness-of-fit test on the basis of the multiplier 
bootstrap we show in Table [6] the simulated rejection probabilities of the pdm-bootstrap test 



(4.9) 



where qf^^^ denotes the (1 — a) quantile of the bootstrap distribution. For the null hypothesis we 
considered as the parametric class the family of Clayton tail copulas. In particular we investigated three 
scenarios corresponding to a coefficient of tail dependence Ai(l, 1) varying in {0.25, 0.5, 0.75}. Under the 
alternative we consider three models: 

(1) A convex combination between the independence copula n(u) = U1U2 and a Clayton copula [with 
convex parameter 1/3], such that the tail copula is given by A2,(x) = 1/3 (x^^ + ^2^)^^^^. The 
parameter 9 is chosen such that Al = Ai(l, 1) = 1/3 x 2~^/^ varies in the set {1/12,2/12,3/12}. 



(2) The asymmetric negative logistic model [see Joe (1990)], defined by 

-i/e 



A 



Xl-t,t) = {{M^-t))-' + {i,2t)-'\ 



t G [0,1], 



with parameters ipi = 2/3, ip2 = 1 and 9 G (0, 00) chosen such that A^, = Aj;,(l, 1) varies in the set 
{0.2,0.4,0.6}. 



(3) The mixed model [see Tawn (1988)], given by 

AL{l-t,t)=9t{l-t), te[0,l], 
where the parameter 9 £ [0, 1] is chosen such that A^ = Ai(l, 1) = 9/2 equals 0.1 or 0.3. 
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n 1 '7/1 
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(_/layton(AL 
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J. / \ 

(_/layton(AL 


= 0.75) 


0.091 


0.048 
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0.091 


0.058 


0.024 


Convex(AL 


= 1/12) 


0.095 


0.052 


0.017 


0.386 


0.291 


0.179 


Lyonvex^Ai 


= 2/12) 


0.124 


0.066 


0.029 


0.502 


0.401 


0.253 


Convex(AL 


= 3/12) 




U.2UU 


U.Uoo 


n QQCi 

U.ooU 


U.ozo 


U. (UU 


Aneglog(AL 


= 0.2) 


0.119 


0.071 


0.028 


0.257 


0.185 


0.109 


Aneglog(AL 


= 0.4) 


0.241 


0.174 


0.105 


0.625 


0.534 


0.416 


Aneglog(AL 


= 0.6) 


0.874 


0.833 


0.732 


1.000 


1.000 


1.000 


Mixed (Al = 


0.1) 


0.118 


0.069 


0.022 


0.523 


0.424 


0.268 


Mixed (Al = 


0.3) 


0.148 


0.068 


0.032 


0.405 


0.315 


0.187 



Table 6: Simulated rejection probabilities of the pdm-bootstrap test (4.9) for the hypothesis (4.7). The 
first three lines are models from the null hypothesis, whereas the last eight lines correspond to alternatives. 
The sample size is n = 1000 and B = 500 Bootstrap replications have been performed. \l denotes the 
lower tail dependence coefficient. 



The results are based on 1000 simulation runs, while the sample size is n = 1000 and two cases k = 
50, 200 are investigated for the choice of the parameter k. For each scenario the critical values have been 
calculated hy B = 500 bootstrap replications with ZY({0, 2})-multipliers. We observe a reasonable power 
and approximation of the nominal level in most cases. Under the null hypothesis the test is conservative 
and this effect is increasing with the level of tail dependence. For the mixed model with /c = 50 the power 
of the test is close to the nominal level. This observation can be explained by the fact that for A^ = 0.5 
[which corresponds to the case 6 = 1] the model is exactly the same as the Clayton model with parameter 
1, i.e. we get close to the null hypothesis with increasing tail dependence. Finally, we note that a larger 
choice of the parameter k results in substantial better power properties, while we do not observe notable 
differences in the quality of the approximation of the nominal level. Again, this may be explained by the 
fact that bias terms in GOFn cancel out when calculating the difference Hn = Vk — Al{-;9^^^)^. 
Therefore, we propose to use rather large values for k in applications of the goodness-of-fit test. 
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A Proofs 



A.l Proof of Theorem O 



Let ;Soo(IR+) denote the set of functions / : M+ — )• M (where M+ = [0,oo)) that are uniformly bounded on 
compact sets (equipped with the topology of uniform convergence on compact sets) and define Bi^{M.^) 
as the subset of all non-decreasing functions / : M-|_ — t- R^- which satisfy /(0+) = and for which 
sup ran/ < oo implies that there exists a xq with /(xq) = sup ran/. The latter condition implies that 
the adjusted generalized inverse function defined by 



sup{x G M+ I fix) = 0}, 

inf{x e M+ I fix) > z}, 

inf{x G M+ I fix) = sup ran/}, 



stays in Ba 



for every / e (M+). Further set 
B'JiRl) := {7 e B^iRl) I 7(.,oo) G ^^"(II 



z = 

< z < sup ran/ 
z > sup ran / 



nioo,-)eB'J>iR+)} 



and now define a map <I> : B^^i 



7 



$(7) 



fioo(t^) by 

7(7-(x,oo),7-(oo,y)) 
7(7" (x, 00), 00) 
7(oo,7~(oo,y)) 



if x,y ^ 00 
if y = 00 
if x = 00, 



see also Schmidt and Stadtmiiller (2006). Observing that A/, G Bi^iM^j^_) and that the adjusted generalized 
inverse of A£,(a;,oo) is given by ^FiiF~^ikx/n)), one can conclude that $(Ai) = Al and <^*(A2,) = A^, 



almost surely) and the proof of Theorem 2.2 follows from the functional delta method (Theorem 3.9.4 



m 



Van der Vaart and Wellner (1996)) and the following Lemma, which is an extension of the result in 



the proof of Theorem 5 in Schmidt and Stadtmiiller ( 2006 ) 



Lemma A.l. Let be a lower tail copula whose partial derivatives satisfy the first order property (2.13) 
for p = 1,2. Then ^ is Hadamard-difjerentiable at Al tangentially to the set 

C°(]R^) = {7 G Soo(K+) I 7 continuous with 7(-,0) = 7(0, •) = O} . 

Its derivative at Al in j € C'^(M^) is given by 

^Ai(7)(x) = 7(x) - 5iAi(x)7(xi,oo) - 92Al(x)7(oo, X2) (A.l) 

where dpA^, p = 1,2 is defined as on the set {x G | Xp G {0, 00}}. 

Proof. Decompose ^ = o ^2 ° ^1 where 

$1 ■.B'J'iRl) ^ B'J>iRl) X B'J>iR+) X ^^°(M+) 



7 

$2 -.B, 



1,0 ( 

00 V 



(7, 7(-, 00), 7(00, •)) 
ID X S^°(M+) X S^°(M+) 



B'JiRi) X Bi^^'-iR^) X Bi^''-iR+) 



?7,0 



7,0, 

00 



iiJ,g) ' — ^ (7,/ ,9 )) 
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i{f{x),g{y)) , if x,y/oo 
{lJ,g)' — ^ S 7(/(2;),oo) , if y = oo 
_7(oo,5'(y)) , if x = oo, 

where ;S(^'^'~(R+) denotes the set of all adjusted generalized inverse functions /~ with / G B^oo(R+)- Now 
<I>i is Hadamard-differentiable at A^, tangentially to C'^(M^) since it is linear and continuous. The second 
map <I>2 is Hadamard-differentiable at (A/,, ids^, id]R_^) tangentially to C'^(M^) x C''(M+) x C''(M+) where 
C'^(M+) consists of all continuous functions / on M+ with /(O) = and its derivative at (A^, idjR^, id]R_^) 
in (7,/, 5) is given by $2 (a^^ idg idg )(')'' fid) — (7) ~fi ~9)- The proof follows along similar as the one 
of Theorem 5 in [Schmidt and Stadtmiiller ( [2006 ), p. 321 and is therefor omitted, we just note that 



1+ +tnfn) (x) > for all X > is implied by the additional assumption of continuity in for functions 
m the set ^^'°(M+). Some more efforts are necessary to show that $3 is Hadamard-differentiable at 
(AL,idK_^,idiR^) tangentially to C°(l^) x C°(M+) x C°(M+) with derivative 

^3,(AL,id.+ ,idM^)(7> />5)(x) = 7(x) + aiAL(x)/(xi) + 92Al(x)5(x2). 

To see this let t„ ^ 0, (7^, /„, gn) G Soo(t^) x 600 x Boo{R+) with (7^, /„, g^) ^ (7, /, g) e CO(l^) x 
CO(M+) X CO(M+) such that (A^ + t„7„,idM+ +tn/n,idK+ +tn9n) G Bh\Rl) X S^°'"(M+) X S4°'-(M+). 
Now $3 is linear in its first argument and we introduce the decomposition 

t~'^ {$3(Al + *n7n,idlR+ +tn/n,idlR+ "htnS'n) " ^3 (Al, id]R_^ , idlR^ ) } = L„i + L„2, 

where 

Lnl = {^3(AL,idK_^ +tn/n,idR+ +tngn) " ^'3 ( Al , idR+ , idjR^ ) } 
Ln2 = $3(7n,idR+ "Ftn/n, idR^ +tngn)- 

By the definition of d it suffices to show uniform convergence on sets T of the form T = [0, Mi] x {cxd} U 
{00} X [0,M2] U [0,M3]2, where Mi,M2,M3 £ N. Since T C is compact {fn,gn) converges uniformly 
and 7 is uniformly continuous; hence Ln2 uniformly converges to 7. 

Considering L^i we split the investigation into six different cases. First, let x E (0,M3]2. A series 
expansion at x yields 

Lnl = 9iAl(x)/„(3;i) + d2AL{x)gn{x2) + r„(x), 

where the error term r„ can be written as 

r„(x) = (aiAi(y) - 5iAi(x))/„(xi) + (52AL(y) - 52Al(x))5„(x2) 

with some intermediate point y = y{n) between x and (xi -|- tnfn{xi), X2 + tnfnix2))- The dominating 
term converges uniformly to 9iAj;,(x)/(xi) -|-92AL(x)(5r(x2), hence it remains to show that r„(x) converges 
to uniformly in x. For a given e > uniform convergence of fn and uniform continuity of / on [0, M3] 
as well as the fact that /(O) = allows to choose a 6 > such that |/n(xi)| < e for all xi < 5. Since 
partial derivatives of tail copulas are bounded by 1, the first term of r'„(x) is uniformly small for xi < 6. 
On the quadrangle [5, M3] x (0, M3] the partial derivative diA^ is uniformly continuous which yields the 
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desired convergence under consideration of y(n) — )• x and boundedness of /. The same arguments apply 
for the second derivative and the case x E (0,M3]2 is finished. 

Now consider the case x G (0, M3] x {0}. By Lipschitz-continuity of Al on we get 

\Lnl{xi,0)\ = t~^|AL(xi +tnfnixi),tngn{0))\ 

= t~^\AL{xi + t„/n(xi),i„5rn(0)) - Al(xi + t„/„(xi),0)| 
< bn(0)| ^5(0) =0. 

Since diAL{xi,0)f{xi) + f?2Ai(xi, O)g(O) = this yields the assertion. For the cases x = (0,0)"^ and 
X G {0} X (0, M3] the arguments are similar and we proceed with x G [0, Mi] x {00} (and analogously 
X e {00} X [0,M2]) 

Lnlixi,(X>) = t~^{AL{xi + t„/„(xi),00) - Al{xi,(X>)) = fn{xi) f{xi). 

By 9iAi(xi,c«) = 1 and 32Al(xi,oo) = this yields the assertion. To conclude, $3 is Hadamard- 
differentiable as asserted. 

An application of the chain rule (see Lemma 3.9.3 in Van der Vaart and Wellner ( 1996| )) completes the 
proof of the Lemma. □ 



A.2 Proof of Theorem ISTTl 



Due to Lemma C.2 in the Appendix B [which is an analogue of Theorem 1.6.1 in Van der Vaart and 



Wellner ( 1996 ) for the case of conditional weak convergence] the proof of conditional weak convergence of 
in 5oo(M^) can be given for each l°°{Ti) separately. To this end we note that every Tj can be written 
in the form T = [0, Mi] x {00} U {00} x [0, M2] U [0, M^f , where Mi, M2, M3 E N, and show conditional 



weak convergence in 1°°{T). Recalling the notation of /n,x(Ui) in (3.3) we can express as 



<(x) = '^VkiAi - A.) = ^ (I - 1) /„ 



and the assertion now follows by an application of Theorem 11.23 in Kosorok (2008). For this purpose 



we show that the assumptions for this result are satisfied. Let = {/n,x : x G T} be a class of functions 
changing with n and denote by 



Fn(u) 



{ui < kM/n or U2 < kM/n} 



M = Ml V M2 V M3 a corresponding sequence of envelopes of J^n- We have to prove that 
(i) {Fn-,Fn) satisfies the bounded uniform entropy integral condition 



limsupsup / J log N{e\\Fn\\Q ^2 , Fn 1 L2{Q)) de < 00, 

n-^oo Q Jo 



(A.2) 



where for each n the supremum ranges over all probability measures Q with finite support and 

\\Fn\\Q,2 = {! Fn{x)^ dQ{x)y^^ > 0. 



(ii) The limit H{x,y) = lim„_j.oo IE[Q;n(x)an(y)] exists for every x and y in T . 
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{in) limsup„_^ooIEi^«(Ui) < oo 

(iv) lim„,^oo IEF2(Ui)I{F„(Ui) > e^} = for all e > 0. 

(v) lim„^oo Pn(x,y) = /o(x,y) for all x,y G R^, where 

p„(x,y) = (E(/„,,(C/i))-/„,y([7i)) 



2^l/2 



(A.3) 



Furthermore, for all sequences (x„)„, (y„)„ in T the convergence Pn(xn,yn) — ^ holds, provided 

/0(Xn,yn) 0. 

(m) The sequence J-n of classes is almost measurable Suslin (AMS), i.e. for all n > 1 there exists a 
Suslin topological space T„ C T with Borel sets Bn such that 

(a) P*(sup,eTinfyeT„ |/n,x(Ul) - /n,y(Ui)| > 0) = 0, 

(b) fn,- '■ [0, 1]^ X T„ — M is S|[o i]2 x i3„-measurable for z = 1, . . . , n. 

In order to prove the bounded uniform entropy integral condition (i) we decompose Tn = U?=i -^n^ with 
-F^'^ ={/£,xGr} and 

f!^:l(U^) = J\^{Uii < kxi/n} 1{X2 = oo}, f'il{Vi) = j|l{f^i2 < kx2/n} l{xi = oo}, 



/n,x(Ui 



^{Uii < kxi/n,Ui2 < kx2/n} I{xi < oo,X2 < oo}. 



The corresponding envelopes of the classes J-n^ are given by 



F^,'\V^) = < kM/n), f(2)(U,) = ^I{U, 



a < kM/n) 



I{Uii < kM/n, Ui2 < kM/n), 



so that F„(Uj) = max^^-^{Fn*^(Ui)}. If we prove that the sequences {Fn' ,Fn') satisfy the bounded 



uniform integral entropy condition given in (A.2), then the condition holds also for {Fn^Fn) by Lemma 
in the appendix and thus the assertion in (i) is proved. We only consider the (hardest) case of Fn • 
Note that = {/„,x, X G [0, Ms]^} = g^^^ ■ Q^^\ where 

/«,x = {n/k)^l'^\{\J,x < kxi/n,Ui2 < kx2/n}, 
Gk'^ = {9n,t = {n/k)'/H{Uij < kt/n] I t e [0,M3]} 

for j = 1,2. Since the functions gn,t are increasing in t the G){' are VC-classes with VC-index 2. Thus 



(A.2). Proposition 11.22 in 



by L emma 11.21 in Kosorok J2008) both c lasses satisfy the bounded uniform integral entropy condition 



Kosorok 



(3) 

( 2008 ) shows that Fn has the same property and by the discussion 



at the beginning of this paragraph (i) is satisfied. 

For the proof of (ii) note that E[a„(x)an(y)] = n/k - C{^)C{^)^ , which converges to 

Ai(x A y) =: H{^, y), since r^C{f)C{^) ^ 0. 
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Regarding (iii) and (iv) we note that EF„(Ui)^ = 2M — ^C{Mk/n, Mk/n), which converges to 2M 
Al{M,M). Further, 

EF2(Ui)I{F„(Ui) > eV^} = [ F^iVi) dF 

J{F„{lJi)>eV^} 
Tl 1 

< -F(-I{Un < kM/n or U12 < kM/n} > e) = 
k k 

for sufficiently large n, such that k > 1/e. For (v) we note that 



Pn(x,y) = (E(/„,x(Ui)-/.,y(Ui)) 



(C(xA;/n) - 2C((x A y)k/n) + C{yk/n)f^ 



(Ai(x) - 2Ai(x A y) + AiCy))^/' =: p(x, y). 



Due to Theorem 1 in [Schmidt and Stadtmuller| ( [2006 ) we have locally uniform convergence in the latter 
expression, which yields the second condition stated in (v). 



For the proof of condition (vi) we use Lemma 11.15 and the discussion on page 224 in Kosorok (2008) 



and show separability of i.e. for every n > 1 there exists a countable subset <ZT such that 

sup inf |/„,y(Ui) - /„,x(Ui)| > 0) = 0. 

Choose T„ = (Q n [0,Mi] x {00}) U ({00} x Q n [0,M2]) U (Q^ n [0,M3]2), then we have (note that 
the functions /„^x are built by indicators) that for every uj and every x G T there is an y G T„ with 



|/n,x(Ui(a;)) — /„^y(Ui(a;))| = 0. This yields the assertion and thus the proof of Lemma 2.1 is finished. □ 
Remark A. 2. Observing that 



a™ X 



1 " 

^ V(/„.x(U,)-E/„,x(Ui)) 

i=l 



and that in Section A. 2 we showed the sufficient conditions for an application of Theorem 11.20 in Kosorok 



(2008), we obtain an alternative proof of Lemma 2.1 



A.3 Proof of Theorem lOl 



For technical reasons we give a proof of Theorem 3.4 in advance of Theorem 3.2 and 3.3, The proof 



is essentially a consequence of a bootstrap version of the functional delta method, see Theorem 12.1 in 



Kosorok (2008). Since this result only holds for Banach space valued stochastic processes some adjust- 
ments have to be made. Note that the space i3oo(]R+) is a complete topological vector space with a metric 
d and some care is necessary whenever technical results depending on the norm are used. 
Due to Lemma 12.11 and Theorem 13.11 we have 

r 

in Soo(M^). Observing that the generalized inverses of Al{x, 00) and A|^(x, 00) are (P-almost surely) given 
by '^Fi{F~^{kx/n)) and ^Fi{F^'[ (kx/n)), respe ctive ly, one can conclude that <I>(Al) = Al,^{Al) = Al 

^ is Hadamard-differentiable on i3^(M^) at 70 = A^, 



k{AL - Al) - , Vk^{Al- Al) 



and $(Ai) = A^'^ (P-almost surely) . By Lemma 



A.l 
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tangentially to C°(M^) C Soo(M^). Therefore it remains to argue why Theorem 12.1 in Kosorok pOOS^ 
can be apphed in the present context. 



A careful inspection of the proof of Theorem 12.1 in Kosorok (2008) shows that properties going beyond 
our specific assumptions (i.e. the complete topological vector space (i3oo(K+), d)) are used only three 



times. First of all the mapping needs to be extended to the whole space Bo 



which is possible 



using equation (A.l ) as the defining identity. Secondly, the proof of Theorem 12.1 in Kosorok (2008) uses 



the usual functional delta method as stated in Theorem 2.8 in the same reference, but this result can 



be replaced by Theorem 3.9.4 in Van der Vaart and Wellner (1996), which provides a functional delta 



method holding in general metrizable topological vector spaces. Finally, the proof of Theorem 12.1 in 



Kosorok (2008) makes use of a bootstrap continuous mapping theorem, see Theorem 10.8 in Kosorok 



(2008), which would yield that 



In our specific context this statement follows immediately from the Lipschitz continuity of the derivative 
<I>^^ and an application of Lemma C.3 in Appendix [Cj □ 



A.4 Proof of Theorem |M1 



Consider the mapping ^ : Bjf{M.%) x Bj;{RX) — > Boo{^i) defined by ^ 
$2 are defined in the proof of Lemma A.l| and ^'i is given by 



$3 o $2 o ^1, where $3 and 



I fit 



B. 



1,0 ( 



X ^^0(M+) X Bi^'{R+) 



(/3,7)^(/3,7(-,oo),7(oo,-)). 



Note that we obtain the representations ^'(Al,Al) = A^, ^(Al,Al) = A^ and ^(A]^,Al) = A 
almost surely). Clearly, ^1 is Hadamard-differentiable at {Al,Al) since it is linear and continuous. $2 
and $3 are Hadamard-differentiable tangentially to suitable subspaces as well, see the proof of Lemma 



A.l By an application of the chain rule, see Lemma 3.9.3 in Van der Vaart and Wellner (1996), we can 



conclude that ^ is Hadamard-differentiable (A^, A^) tangentially to C"(M^) x C"(M^ 
^(Al At) ('^' ^) = - ^1^^- i^hi^i > 00) - 92 Al (x)7(oo, X2). 



with derivative 



(A.4) 



Note that, unlike in the previous proof, we do not have weak convergence (resp. weak conditional con- 
vergence) of \^ (^{Al,Al) — {Al,Al)^ and ^Vk ^(A^, A^) — {Al,Al)^ towards the same limiting field, 
which would be necessary for an application of the functional delta method for the bootstrap [see for 



example Theorem 12.1 in Kosorok (2008)]. Nevertheless, we can mimic certain steps in the proof of this 
theorem to conclude the result. To be precise, note that we obtain by analogous arguments as on page 



236 of Kosorok (2008) that 




■.-1(1 



G-2 



unconditionally, where Gi and G2 denote independent copies of G^^ and c = fir ^ . Hadamard-differentiability 
of the mapping (/3,7) 1— >• (^'(/3, 7), ^'(7, 7), (/3,7), (7,7)) and the usual functional delta method [Theorem 
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3.9.4 in Van der Vaart and Wellner (1996)] yields 



/^(Ai,Ai)-^(Ai,AL)\ 
^(Ai,Ai)-V(AL,AL) 

(Ai,Ai)-(Ai,Az.) 
V (Al,Al) - (Al,Al) J 



/^(A,.A,)(c"'Gi+G2,G2)\ 
%„A,)(G2,G2) 

(c-iGi+G2,G2) 



Observing that 



(Ai,AL) 



is linear we can conclude that 



cyk 



*(Ai,Ai)-^(AL,AL) 
(Ai,Ai)-(Ai,Ai) 



(A.,A.)(^1'0) 

(Gi,0) 



(Gi,0) 



Continuity of the map (a, /3, 7) 1— >• (i(a, /3) yields 

d(cV^ (^(Ai,Ai)-M/(Ai,Ai)) ,cVfc (Ai-Ai)) ^0 

in outer probability and thus by boundedness of the metric d also in outer expectation. Since c\/k(K^^ 



Kl) ^ Gi we obtain the assertion by Lemma 



A. 5 Proof of Theorem 



C.4 



□ 



Let r be a set of the form T = [0,Mi] x {00} U {00} x [0,M2] U [OjMa]^, see also the beginning of the 
proof of Lemma 2.1 We start the proof with an assertion regarding consistency of dpKi and claim that 
for any 5 G (0, 1) 



sup 



apAi(x)-5pAL(x] 







(A.5) 



in outer probability. For a proof of (A.5) split T into three subsets as indicated by its definition and then 



proceed similar as in the proof of Lemma 4.1 in Segers (2011). The details are omitted. Regarding the 
assertion of the Theorem we set 

af^ (x) = (x) - ai Ai (x)/3„ (xi , 00) - 52 Ai (x) /3„ (00 , X2 ) . 

Under consideration of Lemma 



C.4 



it suffices to prove that (i(a„ 0™"'") converges to in outer proba- 



pdm -pdm 

n 1 Oln 

bility. By the definition of d we have to show uniform convergence on the set T. Since laff"^ -pdm 
Dni + Dn2, where 



D 



nl 



\M;oo)\, 



D 



n2 



d2AL - d2AL 



l/3n(00,.)| 



we can consider both summands D^p separately and deal with Dni exemplarily. First consider the case 
X e [0, Ms]^, then for arbitrary e > and 5 £ (0,1) 



sup -Dni(x) > e < 

.x6[0,A/3]2 / 



sup Dni{^) >e/2\ + : 
,xe[o,A/3]2,xi><s / 



sup Dni{^) > e/2 

. xe[0,A/3]2,xi<(5 y 



(A.6) 
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Since diKi is uniformly consistent on {x G [0, Ma]^ \ xi> 5} and since fin is asymptotically tight in 1°°{T) 



[Pn converges unconditionally by the results in Chapter 10 of Kosorok (2008)] the first probability on the 
right-hand side converges to zero. 

Regarding the second summand note that F~^{kx/n) = (where [x] = min{A; € Z | A; > x}) so 

that 



sup 

xe[o,M3]2 



5iAl(x) 



< 



sup 

xe[0,M3]2,3;i>/i 



\k{x^ + h)^-\k{x,-h)^ M3 

2h - 2kh - 



for sufficiently large n. Hence the right-hand side of equation (A.6) is bounded by 



F* sup |/3„(x)| > e/4 , 

yxe[0,A/3]2,xi<5 J 

eventually. As ~^ (unconditionally) the lim sup of this outer probability is bounded by 



sup I 

,xe[0,M3]2,a'i<<5 



^A,(x)| >e/4 



Since G^^ has continuous trajectories and Ga^(0,X2) = (almost surely) this probability can be made 
arbitrary small by choosing 5 sufficiently small. The case x G [0, M^]"^ is finished. For x € [0, Mi] x {cxd} 
the arguments are similar, while for x G {oo} x [0, M2] we have Dni = and nothing has to be shown. 
To conclude, sup^g^ -Dni (x) converges to zero in outer probability and because the term sup^^rp can 
be treated similarly the proof is finished. □ 



B Partial derivatives of tail copulas 

Proposition B.l. The first partial derivative of a (lower or upper) tail copula A satisfies 



diA{0,x) 



limt_^oo A(l, t) if X £ (0, 00) 
ifx = 0. 



As a consequence, the only tail copula that admits for continuous partial derivatives in the origin is the 
tail copula corresponding to tail independence, i.e. A = 0, for either the lower or the upper tail. 



Proof. By groundedness and homogeneity of A, see Theorem 1 in Schmidt and Stadtmiiller (2006), we 
have 

A(/i,x) - A(0,x) 



9iA(0, x) = lim 



h 



lim k{l,x/h) = lim A(l,t) 



for all X G (0, 00). Similarly, OiA(0,0) = 0. The addendum follows by Theorem 1 {iv) in Schmidt and 



Stadtmiiller (2006). 



□ 



As an example, note that for the Clayton copula given in (3.11) we obtain 9iA2,(0,x) = 1 for all ^ > 0. 
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C Auxiliary results 



In the last section we present several technical details. We omit the proofs of the assertions and refer the 



reader to the thesis Biicher (2011). 



Lemma C.l. Suppose Qn and Tin are sequences of measurable functions with envelopes Gn and Hn, so 



that {Gn,Gn) and (T-LnjHn) satisfy the bounded uniform integral entopry condition as stated in (A. 2} 



Then the bounded uniform entropy integral condition (A. 2} holds also for Tn = Gn^ T~in, with envelopes 

Fn = Hn ■ 

Lemma C.2. Suppose Gn = G'„(Xi, . . . X„, ^i, . . . ^n) is some statistic taking values in i3oo(M+). Then 



a conditional version of Theorem 1.6.1 in Van der Vaart and Wellner (1996) holds , namely Gn G in 



;Soo(M+) is equivalent to Gn G in l°°{Ti) for every i G N. 

Lemma C.3. Suppose that g : Di — > O2 is a Lip schitz- continuous map between metrized topological 
vector spaces. If Gn = Gn{Xi, . . . . . . G in Di, where G is tight, then g{Gn) g{G) in 

1D)2. 

Lemma C.4. Let Yn = • • . , 6, • • • , Cn) and Zn = Zn{Xi, ... , X„, ^i, . . . , Cn) be two (bootstrap) 

statistics in a metric space {D,d), depending on the data Xi, . . . ,Xn and on some multipliers ^]^, . . . , i^^. 

F P P 

IfYn-^Y in where Y is tight, and d(Yn, Zn) — >■ 0, then also Zn Y in D. 
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